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Interaction blockade occurs when strong interactions in a confined few-body system prevent a 
particle from occupying an otherwise accessible quantum state. Blockade phenomena reveal the 
underlying granular nature of quantum systems and allow the detection and manipulation of the 
constituent particles, whether they are electrons [l], spins [2], atoms [3HS]? oi" photons 0. The di- 
verse applications range from single-electron transistors based on electronic Coulomb blockade^ to 
quantum logic gates in Rydberg atoms 8 , 9 . We have observed a new kind of interaction blockade 
in transferring ultracold atoms between orbitals in an optical lattice. In this system, atoms on 
the same lattice site undergo coherent collisions described by a contact interaction whose strength 
depends strongly on the orbital wavefunctions of the atoms. We induce coherent orbital excita- 
tions by modulating the lattice depth and observe a staircase-type excitation behavior as we cross 
the interaction-split resonances by tuning the modulation frequency. As an application of orbital 
excitation blockade (OEB), we demonstrate a novel algorithmic route for cooling quantum gases. 
Our realization of algorithmic cooling[10 , 11 utilizes a sequence of reversible OEB-based quantum 
operations that isolate the entropy in one part of the system, followed by an irreversible step that 
removes the entropy from the gas. This work opens the door to cooling quantum gases down to 
ultralow entropies, with implications for developing a microscopic understanding of strongly corre- 
lated electron systems that can be simulated in optical lattices 12", T3^. In addition, the close analogy 
between OEB and dipole blockade in Rydberg atoms provides a roadmap for the implementation of 
two-qubit gates ^2 in a quantum computing architecture with natural scalability. 



An ultracold gas of bosonic atoms in the ground band 
of an optical lattice is described by the Bose-Hubbard 
model [15], in which atoms can tunnel between neighbor- 
ing sites and interact via an onsite repulsive contact in- 
teraction. In a deep lattice where the interactions domi- 
nate, the ground state of the system is a Mott insulator 
with a fixed atom number per site that is locally con- 
stant over a region of the insulator [16]. The energy per 
site in the absence of tunneling is ^Uggn{n — 1), where 
Ugg is the interaction energy for two atoms in the ground 
lattice orbital state and n is the atom number on the 
site. The Mott state exhibits a transport blockade phe- 
nomenon in which the presence of an atom on a site 
energetically prevents tunneling of a neighboring atom 
onto that site even in the presence of a small bias be- 
tween the sites. The transport is blocked unless the bias 
makes up for the interaction cost, making it possible, 
for example, to count atoms tunneling across double- 
wells in a superlatticepl. In this work, we explore an 
excitation blockade phenomenon that does not involve 
transport in the lattice. The excitation transfers local- 
ized atoms between different orbitals on the same site 
through modulation of the lattice depth at a frequency 
close to a vibrational resonance. Physics in higher optical 
lattice orbitals has been the focus of much recent exper- 
imental work including the study of dynamics in higher 
orbitals pT]. multi-orbital corrections to the interaction 
energy [T8^, and unconventional forms of superfluidity in- 
volving higher orbitals[19, 20 . 

The OEB mechanism can be understood in the sim- 
plest scenario for two atoms in a single site of a deep 
three-dimensional lattice, in which the vibrational fre- 




FIG. 1. Orbital excitation blockade mechanism in an optical 
lattice. A single atom on a site is excited to a higher orbital 
by resonantly modulating the lattice depth. For two atoms on 
the same site, interactions lead to an orbital-dependent energy 
shift. Modulation at the appropriate frequency excites one of 
the atoms to the higher orbital, but is off-resonant for exciting 
the second with a blockade energy 6. 



quencies in all three directions are taken to be different 
to avoid degeneracies. The lattice depth along the z- 
direction is modulated weakly, which in the presence of 
anharmonicity of the lattice potential drives atoms be- 
tween the ground orbital and a single, specific excited z- 
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FIG. 2. Time, frequency and site-resolved coherent transfer 
of atoms in a Mott insulator between orbit als. a, Excitations 
transferring a single atom in the n = 1 (orange) or n = 2 
shell (blue) from the ground to the second excited orbital 
are spectroscopically resolved in a two-shell Mott-insulator. 
b, Rabi oscillations between the two orbit als are observed by 
driving at the resonant frequencies for atoms in the n — 1 shell 
and c, n = 2 shell of a Mott insulator. Bose-enhancement 
leads to faster oscillations in the n — 2 shell. When the 
atom number is reduced to obtain one atom per site in the 
region previously containing two atoms, the interaction shift 
suppresses oscillations (black) . All error bars are one standard 
error of the mean, (i-iii) Site-resolved snapshots of the Mott 
insulator are shown at different points in the Rabi cycles. 



orbital, subject to a selection rule that only allows cou- 
pling to orbitals of the same symmetry. For a single 
atom, excitation to the mth orbital requires modulation 
at a frequency ujz^o^m which is approximately mujz^o^i 
ignoring the anharmonicity of the onsite potential. With 
more than one atom on a site, the interaction introduces 
an orbital-dependent shift of the energy levels as shown 
in Fig. [l] In general, the interaction shifts /7gg, /7ge, and 
Uee are all different and the differences are a significant 
fraction of /7gg, where g{e) denotes atoms in the ground 
(excited) orbital. If the coupling strength due to the 
modulation is small compared to these differences and the 
modulation frequency is tuned to uJz^o^m-\-{Uge — Ugg)/h^ 
only a single atom is transferred to the higher orbital 
and the transfer of a second atom is off-resonant. In this 
sense, the first excitation blocks the creation of a second 
excitation. 

The experimental system has been described in pre- 
vious work [2 1. A two-dimensional Bose-Einstein con- 
densate of rubidium atoms resides in a single plane 
of a one-dimensional optical lattice, henceforth referred 
to as the axial lattice, with a vibrational frequency of 
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FIG. 3. Algorithmic cooling in an optical lattice, a, Landau- 
Zener chirp for transferring entropy from the ground to the 
fourth band. The lattice modulation frequency is swept across 
the transition resonances from left to right. The interaction 
shifts Auz for excitation of one of n atoms to the fourth or- 
bital relative to the excitation frequency for a single atom on 
a site, are shown for different orbital occupations. Excitation 
processes in the same column happen at almost the same fre- 
quency to within 30 Hz (see Supplementary Table [l|. b, A 
state with random occupation in a deep lattice is far from 
the Mott insulating ground state. Sequential filtering opera- 
tions followed by reduction of the confinement prepares the 
ground state, which can be adiabatically converted to a ther- 
malized superfluid in a shallow lattice. Red (blue) spheres 
denote atoms in the ground (excited) band. 



^z,o^i = 27r X 5.90(2) kHz. The z-axis is perpendicular 
to the plane and points along the direction of gravity. 
In addition, we introduce a lattice in the plane with a 
spacing of a = 680 nm and a depth of AbEr (trap fre- 
quency of 17 kHz), where Er = h'^/Sma^ is the recoil 
energy of the effective lattice wavelength, with m the 
mass of ^^Rb. The resulting Mott insulator is at the fo- 
cus of a high resolution optical imaging system capable 
of detecting atoms on individual lattice sites through flu- 
orescence imaging. Light-assisted collisions at the start 
of the imaging process reduce the occupation of a site to 
its odd-even parity [2 1. 

We start by demonstrating coherent driving of atoms 
in a Mott insulator between two orbitals. In the presence 
of a harmonic trap, the atoms in a 2D Mott insulator 
are arranged in concentric rings of fixed atom number 
per site, known as shells, with the largest occupation at 
the center [21]. We prepare a Mott insulator with two 
shells and modulate the axial lattice depth by ±1.1(1)% 
at a frequency chosen to transfer atoms from the ground 
state to the second excited orbital. A modulation fre- 
quency corresponding to exactly CjOz^o^2 is resonant for 
atoms in the outer shell with one atom per site (n = 1). 
Excitation to the fourth excited orbital is suppressed be- 
cause of an energy shift of h x 1200(80) Hz owing to the 
anharmonicity of the onsite potential. Rabi oscillations 
between the states \g) and |e) are detected by lowering 
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the axial lattice depth at the end of the modulation so 
that the excited orbital state becomes unbound and any 
population in it escapes along the z-axis due to gravity. 
The Rabi oscillations in that shell, shown in Fig.[2]3, have 
a frequency VL = 2ttx 23.3(2) Hz. 

The OEB is demonstrated in the inner n = 2 shell by 
modulating at a frequency of oj z, 0^2 {Uge — Ugg)/fi. 
For our parameters, /7gg, /7ge, and Uee are 
h X 480(30), 360(20), 310(20) Hz respectively. The 
Rabi oscillations between \g^g) and l/>/2(|e,^) + |^,e)) 
are detected as an oscillation of the parity between even 
and odd after ejecting the atom in the excited orbital 
and are shown in Fig. For the same modulation 
amplitude as before, the oscillations are expected to 
occur ^/2 times faster than the resonant oscillation in 
the n = 1 shell owing to Bose-enhancement[5j. We 
indeed observe a frequency ratio of 1.42(1) between the 
oscillation frequencies. A full frequency spectrum in 
the two shells is shown in Fig. [2^ and the frequency 
separation of the two resonances of 160(10) Hz matches 
well with the theoretically expected value of 165(35) Hz 
when the impact of virtual orbital changing collisions is 
included (see Methods). 

We next employ OEB to demonstrate a new path 
to cooling quantum gases. Evaporative cooling has 
been the workhorse technique for cooling atomic gases 
to nanokelvin temperatures. However, current interest 
in studying the physics of strongly-correlated materials, 
such as high-Tc cuprates, using ultracold gases [12 , 13 has 
spurred research into developing new cooling techniques 
that can reach the requisite picokelvin regime [22l [23] . 
The field of quantum information offers an alternative 
cooling paradigm, wherein a sequence of unitary quan- 
tum gates purifies a subset of the qubits in a sys- 
tem by moving entropy and isolating it in another part 
of the system[10 . One realization of such a cooling 
scheme, heat-bath algorithmic cooling, has been suc- 
cessfully demonstrated with solid-state nuclear magnetic 
resonance qubits [11. We introduce an analogous tech- 
nique for quantum gases where the unitary operations 
are achieved using OEB, building on previous theoretical 
proposals in this direction[24-27 . 

A bosonic quantum gas at a finite temperature T adia- 
batically loaded into the ground band of a optical lattice 
stores its entropy in the form of atom number fiuctua- 
tions in the zero-tunneling limit. Within the local density 
approximation, a lattice site with a local chemical poten- 
tial /i is described by a density matrix p = Xln^^l^)(^l' 
where p^, the probability of having n atoms on the site, 
is given by e-^(^^^^^(^-^)-^^)/Z. Here, p = l/ksT and 
Z is the grand canonical partition function. Cooling to 
zero temperature is achieved by changing the atom num- 
ber distribution on each site to obtain Pn = ^nj/i/c/^g] • 

The crucial ingredient for algorithmic cooling is a uni- 
tary operation that realizes the transformation |n,m) 



|n — l,m + 1) for each n separately, where |n,m) de- 
notes a Fock state with n atoms in the ground band 
and m atoms in an excited band. Resonant lattice mod- 
ulation in the presence of OEB results in a rotation 
gate Rnm{t) = ex.p[i{Qrimt\n-l,m^l){n,m\^c.c.)], 
where ftnm is the transition's Rabi frequency and the re- 
quired transformation is achieved for a modulation time 
t = 7r/2Qnm- Entropy is transferred from the ground 
band to the excited band by performing a sequence of 
TT gates Rn-s^s from s = to s = N — 1 with N cho- 
sen large enough such that pn ^ in the initial state. 
At the end of this sequence, most of the entropy of the 
gas has been transferred to the excited band and can be 
removed from the system by ejecting the atoms in that 
band. The local chemical potential ja is readjusted to 
recover a situation closer to thermal equilibrium by re- 
ducing the harmonic confinement to a new value cjiow so 
that /i < Ugg throughout the gas. At this point, residual 
entropy is stored in the resulting n = 1 Mott insulator 
in the form of holes that are preferably located near the 
edge of the cloud. The gas is allowed to rethermalize 
by lowering the lattice depth to allow tunneling, and the 
final entropy of the thermalized state would be signifi- 
cantly reduced compared to the initial state. 

We start by experimentally demonstrating the algo- 
rithm on a state with known atom number, namely a four 
shell Mott insulator, and reducing the site occupation ev- 
erywhere in the insulator to a single atom per site. To 
increase the blockade energy for this set of experiments, 
we transfer atoms to the fourth axial orbital rather than 
the second. We also replace the rotation gates demon- 
strated in the first part of this work with Landau-Zener 
transitions to improve the fidelity of the algorithm. We 
linearly sweep the modulation frequency from 20.90 kHz 
to 21.65 kHz in 250 ms. The chirp realizes a sequence 
of quantum operations that transfer atoms to the ex- 
cited orbital one at a time, until only one atom remains 
in the ground state in all shells (Fig. |3|i). We probe 
the ground orbital occupancy at different points in the 
frequency chirp by ejecting atoms in the higher orbital 
as before and then performing the parity imaging. The 
parity of the different shells during the chirp is shown in 
Fig. |4^, and an analogue of the typical Coulomb blockade 
staircase is seen in the data. Shell-sensitive manipulation 
of a Mott insulator had been achieved in previous exper- 
iments using a microwave transition between hyperfine 
states, but the lack of a strong blockade allowed transfer 
of only a small fraction of the population to the target 
state [28 . 

Next, we demonstrate cooling by performing the algo- 
rithm on a state that is far from the many-body ground 
state. To prepare such a state, we non-adiabatically load 
a condensate into a deep lattice, projecting the wave- 
function onto a state with Poissonian site occupancy that 
rapidly loses coherence between sites. Using the same op- 
eration sequence as before, we progressively reduce the 



4 




_1j . i L 1 — \ . \ J 0.1 ^ ' \ ' 1^ ' 1 ^ ^^^^^^^^ 

21000 21200 21400 21600 21000 21200 21400 21600 Recooled superfluid 

Modulation chirp endpoint (Hz) Modulation chirp endpoint (Hz) 



FIG. 4. Experimental realization of algorithmic cooling, a, By chirping the modulation towards higher frequencies, atoms 
in a four-shell Mott insulator are sequentially excited one at a time to the fourth orbital. Population in the higher orbital is 
subsequently ejected at the end of the chirp. The average parity in the n = 1 (orange), n = 2 (blue), n = 3 (black) and n = 4 
(green) shell is shown at different points in the chirp, together with single shot images, illustrating the conversion to a three, 
two and finally one shell insulator, b. The same frequency chirp algorithmically cools a state with random occupancy into an 
n = 1 Mott insulator, observed as an enhancement in odd occupancy. All error bars are one standard error of the mean, c, (i) 
An incoherent cloud does not exhibit an interference pattern in a 5 ms time of flight expansion after adiabatically lowering the 
lattice depth, (ii) Cooling converts the incoherent cloud to a Mott insulator in the deep lattice. After adiabatically lowering 
the lattice depth, a superfluid forms and an interference pattern is obtained in the expansion images. 



randomness of the ground band occupancy, preparing a 
single-occupancy Mott insulator (Fig. [sJd). We enhance 
the odd occupancy from 0.45(1) to 0.76(2) (Fig. [4]d), 
demonstrating significant atom number squeezing limited 
by the conversion efficiency of the Landau- Zener transi- 
tions. To complete the algorithm, we readjust the har- 
monic confinement to obtain a state close to the many- 
body ground state. We verify the ground state character 
by ramping back adiabatically to a 5Er lattice in 100 ms 
and releasing the atoms from the lattice. Without cool- 
ing, we obtain a featureless cloud shown in Fig. [4]3(i), 
indicating an absence of the coherence expected in the 
superfluid ground state. With cooling, the Mott insula- 
tor is adiabatically converted to a superfluid, giving rise 
to matter wave interference peaks shown in Fig. |4]3(ii). 

We now discuss the limits on the entropies that can 
be achieved with algorithmic cooling. The conversion 
efficiency to a single-occupancy Mott insulator is tech- 
nically limited in our system by heating due to sponta- 
neous emission during the sweep and to a lesser extent, by 
the efficiency of the Landau-Zener sweep for clouds with 
large average occupancies (see Methods). While we have 
demonstrated cooling of hot clouds, the single-occupancy 
probability we have achieved using algorithmic cooling in 
a two-shell Mott insulator is 0.94(1). This is comparable 
to what had been previously achieved with evaporative 
cooling, corresponding to an average entropy of 0.27/cb 
per particle [21 . Nevertheless, lattice heating can be 
made negligible by using a further-detuned lattice (e.g. 
1064 nm), while shaped pulses can improve the Landau- 
Zener transfer efficiency [25]. More fundamentally, the 



single-shot cooling algorithm we have implemented is lim- 
ited by initial holes in the Mott insulator which cannot 
be corrected. However, repeated iterations of the algo- 
rithm can circumvent this problem and bring the cloud 
to zero entropy with quick convergence [25 . The cycle 
alternates between (a) using OEB to produce a reduced 
entropy n = 1 insulator in a harmonic confinement cjiow 
(demonstrated above) and (b) adiabatically increasing 
the confinement to cjhigh in the presence of tunneling to 
move hot particles to the center of the cloud where they 
can be removed again. Alternatively, the outer edge of 
the cloud containing the holes can be removed using the 
high resolution available in our system ^29^. 

In conclusion, we have observed a new blockade phe- 
nomenon in optical lattices when exciting atoms to higher 
orbit als, analogous to dipole excitation blockade in Ryd- 
berg atoms. The blockade permits deterministic manipu- 
lation of atom number in an optical lattice. We have used 
it to convert a multi-shell Mott insulator into a singly- 
occupied insulator with over a thousand sites, the largest 
quantum register achieved so far in an addressable sys- 
tem. The same technique allows initialization of regis- 
ters in longer wavelength lattices where a Mott insula- 
tor cannot be prepared [30]. OEB also opens a route to 
implementing quantum gates in optical lattices. Single- 
site addressing [29 , possible with our microscope, can 
perform rotations of individual orbital-encoded qubits 
rather than the global rotations demonstrated in this 
work. Controlled-NOT gates can be implemented by 
conditionally moving the control qubit onto the target 
qubit site, and performing an interaction-sensitive ro- 
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tation of the target qubit[TT. Finally, the algorithmic 
cooling technique we have developed could potentially 
achieve the ultralow entropies required for quantum sim- 
ulation [121 E] and computation in optical lattices, and 
establishes a bridge to quantum information for import- 
ing novel ideas for cooling quantum gases. 
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METHODS 



State Preparation 

Our experiments begin with a degenerate 2D Bose gas 
of ^'^Rb atoms prepared in the |F = l,m/ = — 1) state in 
a single layer of a ID optical lattice with spacing 1.5 jam, 
in the focal plane of a high resolution imaging system as 
described in previous work. The atoms are then loaded 
into a 2D optical lattice with spacing 680 nm, which is 
ramped up to a depth of AbEr adiabatically on either 
a single- or many- body timescale, depending upon the 
experiment to be performed. 



Higher Band Removal 

The orbital blockade is observed through the deter- 
ministic removal of atoms in higher bands of the 1.5 jam 
lattice. For removal of atoms from the second band, this 
is achieved by reducing the depth of this lattice to 3.8 
kHz from an initial depth of 35.8 kHz. Gravity pro- 
duces a shift of 3.2 kHz per well, thus inducing second 
band atoms to Zener tunnel away within a few ms. The 
Landau-Zener tunneling rate from the ground state is 
given by Tlz = ^e"^-/^ ^ 12 Hz. Here g is the grav- 
itational acceleration and gc = auj\l^. This effect leads 
to a loss of ground state atoms on the percent level, but 
can be made negligible by using excitations to the fourth 
band. 
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TABLE 1. Frequency shifts in Hz for a transition transferring 
an atom from the ground band to the fourth band, starting 
with n and m atoms in each of these bands. The shifts are 
measured relative to an initial state with one atom in the 
ground band and none in the excited band. 



up to a combinatoric factor Cmn^pq resulting from 
Bose enhancement. For an energy defect of 
^mn^pq ^ ^mn^pq^ this proccss produccs an energy shift 
of Amn -\^mn^pq\'^ / ^mn^pq- For our experiment, the 
dominant band changing collision is|m = 0,n = 2)^ 
\p = l,g = 1) with a Rabi frequency ^02^11 = 27r x 
120 Hz. In this case lattice anharmonicity and interaction 
shifts produce an energy defect of ^02^11 = 27r x 330 Hz, 
resulting in an additional overall shift of the |0,2) state 
of A02 = -27r X 45 Hz. 



Band Dependent Energy Shifts 

Due to its large spacing, the 1.5 /im lattice has a recoil 
energy of only 250 Hz. Consequently its depth is ~ 150 
recoils, and its low-lying eigenstates are, to good approx- 
imation, the Hermite-Gaussians of a harmonic oscillator. 
The interaction energy between particles in bands m and 
n may thus be written in terms of the ground band in- 
teraction energy Uqo as: 



^00(2 — Snm) 



J\^PM\^\Mx)?dx 
I\Mx)\^dx 



where t/Jmi^) is the normalized mth harmonic oscillator 
wavefunction. The total interaction shift for M particles 
in band m, and TV particles in band n is thus: 



M(M-l)^^ N{N-1)^^ ^^^^^^ 



The interactions also produce ofF-resonant band changing 
collisions, with a Rabi frequency: 



mn-h^pq ^00 ^mn-f-)-pg 



/ 7/^^ (x) (^) V'p (^) (^) 



Limits on entropies achievable with algorithmic 
cooling 

In our system, the main limitation on the achievable 
entropies with algorithmic cooling is losses during the 
Landau-Zener chirp. Spontaneous emission after absorp- 
tion of photons from the lattice leads to excitation of 
atoms from the ground state at a rate of 0.1 s~^. These 
atoms are quickly lost due to tunneling in the higher 
bands in the case of the in-plane lattice or in the band 
filtering step in the case of the axial lattice. This leads 
to holes in an n = 1 Mott insulator (3% during the 250 
ms ramp), setting a lower bound on the reachable en- 
tropy after thermalization of 0.15A:b per particle. The 
lattice light is detuned 25 nm to the blue of the atomic 
resonance, and the heating rate can be made negligible by 
increasing the detuning (e.g. using 1064 nm lattice light). 
A more fundamental limit on the single-shot cooling al- 
gorithm demonstrated here is given by initial holes in the 
Mott insulator that cannot be corrected. For our Mott 
insulators, the hole fraction is on the order of 0.5% cor- 
responding to a post-thermalization entropy of ~ O.OG/c^ 
per particle. This limit can be overcome by the iterative 
algorithm described in the text. For Mott insulators with 
large initial atom numbers per site, it is also important 
to take into account the efficiency of the Landau-Zener 
chirp. For excitation to the second band, the measured 
efficiency is 0.94(1). 



